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ABSTRACT 


This paper describes an efficient method for the numerical] 
evaluation by a high speed digital computer of: the integral of an un- 
correlated elliptical Gaussian distribution over a circle centered at the 
mean of the distribution, |'(A,¢c); the integral of a circular Gaussian 
distribution over a circle offset from the mean of the distribution, P(k, d). 


The methods are programmed for both the NORC and the IBM 
7090. They yield 6 decimal digit accuracy with an average computation 
time of 10 milliseconds per case on NORC and 6 milliseconds on the 
7090, 


Two extensive inverse tables are included. One gives the radius 


Kas a function of | and c, and the other gives the radius 4 as a function 
of P and d, 


ii 
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|. INTRODUCTION 


A problem that often arises in applications of statistical theory 
is that of evaluating tne integral over a circle offset from the origin 
with center (4,4) and radius & of an uncorrelated bivariate normal dis- 
tribution with its mean point at the origin and standard deviations o, 
ando, inthe » and y directions respectively, Oxy being a rectangular 
cartesian coordinate system. 


In problems of military operations research, as discussed in de- 
tail in [8], this integral represents the kill probability, or probability 
of damage (depending on whether the.criterion of interest is total de- 
struction or infliction of a stated degree of damage) of a point target, 
by a single weapon (projectile, bomb, etc.), in the case in which the 
distribution of shots is as stated above and the targetis at (4.41, In 
the terminology of [8], ® is the damage radius of the weapon relative 
to the target, and the conditional damage function is the "cookie-cutter" 
function. The assumption is made that the probability of damage (total 
or of preassigned degree) is unity for impact points (1.3) within the 
circle of radius A and center (4,4) and is zero for impact points out- 
side this cirele, 


While it is clear that the adoption of the "cookie-cutter" condi- 
tional damage function is an idealization, it nevertheless often serves 
as avery useful approximation to reality. This idealization has been 
assumed in many studies involving cither the distribution described 
above or one of its special cases as discussed below, and has been used 
in a number of published tables, for example, [2], [3], [4], [6], [9], 
{i1], (12], [13], [14]. Other types of conditional damage functions are 
discussed in [8]. 


This probability can be written in rectangular coordinates as 


hak fea ht wtr- hy? 


ho oh ok &k 1 l x¥ \¥ 
Hitete, Wipes eo gait Sede ie pas dxdy (1) 
O, O, OY a, 2740). 2 Ox ay: 


heh Jee ft? 
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In polar coordinates, equation (1) becomes 


So, K ,2Qr 7 e 
A Rh & 1 1 h+reoos 0 kersin OY l 
Pl-,-,-—,—) = exp { -— ———_—_-_ + [————___ rdrdQ (2) 
o, G,. 0, 9, 270,90, 2 oy o.. j 
oes i “ “a o . 
where 
x-h=rcos0,  yok=rsin), C2rek. Cen s oe (3) 


The integral in this form is, in general, not integrable analstically. 
Its numerical evaluation however can be accomplished efficiently by a 
method of quadrature. The procedure is described in NWL Report 
#1710, [12], and in [15]. 


Two particular cases of the integral given above are specified by 


(i) setting 4 and 4 equal to zero, (See Appendix D, page 1). 

(2) setting «, equal to «,, (See Appendix E, page 1). 
In both cases the resultant integral is amenable to evaluation by com- 
puting terms of a Taylor series or of an asymptotic expansion. The 
purpose of this report is to exhibit these series, in terms of recurrence 
relations, and to describe their application towards an efficient numeri- 
cal calculation of the integral in each of these two special cases. 


Il. BASIC EQUATIONS 


If 4 and 4 are set equal to zero, then equation (2) reduces to 


as K 
RR = 1 1 . 
P{—,—, 0 0 VCK,e) = — exp ——r? [Rad cos O] $ rdrdd (4) 
7, 9% me = 
7 o o 


where a, ~, 
) ; N = 
Osez—¢l Kem (5) 
a, Oy 
le? lec? 
{=:—— Re (6) 
Ie? ae 
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If o, is equal to o,, then equation (2) becomes 


Tee K pan 
i A ho ok i 1 

(ee *) -rannez| | v0} [(d + rome + sn] sar (7) 
T, Ty Fy Fy 20 : 2 


o 


where the symbol ?(k,d) indicates that P is a function of only two var- 
iables, R andd. The definitions of Rk and d? are: 

._ Pp ae 

k = kK [oy d = ———- (8) 
Since o, and o, are equal, the distribution is circular, the circle of 
integration may always be taken as offset along the x axis as indicated 
by equation (7). P(&,d) is known as the circular coverage function. 


Equations (4) and (7) can also be written in the following form: 


ie 
B Ar? 
| exp (-}) 1, (=) rdr (9) 


9 


V(RK,¢) = 


7 fo 


h 
PCR, d) ~ onine 7 | exp(—r? /2) 1, (rd) rdr (10) 


oO 


where I,(x) is known as the modified Bessel function of the first kind of 
zeroth order. ‘This function can be defined by the integral relation 


fn 
il 
I(x) =| exp (v cos 0) dd (11) 
0 


7 


or by its Maclaurin expansion as given in equation (23) below, (see [5]). 


If ¢- 0 in equation (10), the integral reduces to 


Rk 
i exp(-r?/2)rdr, since [,(0) = | (put x = 0 in equation (23) in section V). 
a 


If « - | in equation (9), then 4 =0, B= 1, from equations (6), and the 


eK 
integral becomes J exp(-—r?/2)rdr . These integrals are evaluated 
oO 
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analytically inthe forms [ 1 — exp(-A?/2) ] and [ 1-— exp(-K?/2) ] re- 
spectively. These special cases were pointed out by Germond in [4] 
(for P(A, d)) and Esperti in [2] and Fettis in [3] (for V(K,c)). Another 
special case which can be mentioned here, although it requires the use 
of an equation of Section IV, is that in which k =d in the circular 
coverage function. In this case the V term in equation (18) vanishes, 
and P(R,d) is given analytically in terms of an exponential and a Bessel 
function. This was also pointed out by Germond in [4]. 


These special cases, however, do not require any modifi- 
cations of the computing programs herein described. For the case 
e=1 inthe V(kA,c) function, see the analysis following equation (36) 
in Section V. The case d= 0 is analogous in the P(Rk,d) function, in 
view of relation (74) in Section VI, and can be analyzed similarly, but 
the details are not given inthis report. For the case & =d in the 
P(h,d) function, see the analysis following equation (95) in Section VI. 


If ¢ = 0 inthe V(k,c) function, implying that Oy = co, = 0 (see 
equations (5) ),the distribution is not a bivariate or planar one, and 
the integrals in equations (1), (2), (4) and (9) cannot be evaluated as 
they stand, A and # of equations (6) being undefined. Nevertheless, 
for very small c, the I’(k,¢) function is approximated by a one- 
dimensional probability integral 


> pki 
Erf(K /f'2) = al exp(—y? ) dy (see appendix B). It can be proved 
7 
oO 


rigorously that I'(A,c) tends uniformly to Erf(K/V2)for 0<K <m ase 
tends to zero; that is, for arbitrary «>0, there exists a number 

c,(0 < ¢, < 1), depending only on ¢, such that 6 < Frf(t /¥2) - V(K,0) <e 
for all ¢ such that 0<c¢ <¢, and all positive A. Thus J’(k,0) may be 
consistently defined as Frf(K/¥2). Fettis in equation (I-17) of [3] 
states an equivalent result. Moreover, if |(4,0) is defined in this 
way, it is computed by the regular procedure described in Sections V 
and VI, and hence need not be regarded as a special case from the 
computational point of view. See footnote to equation (50). 


The quantities a! /dkK, and op /dk are used for the inverse 
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table computations which are discussed in Section VII; they are given 
by: 


OV UO«K ‘AK? 

eenars —BK? / ai (12) 
ee exp ( /2) ne ) 

ap R? yd? 

2 tool a ) 1,( Rd) (13) 
OR 2 


Ill. BACKGROUND 


An extensive table of K as a function of V and e had been com- 
puted on NORC, at the request of Dr. H. Weingarten, Special Projects 
Office, Washington, D. C., using the methods advocated herein [10], 
[11] for this report that inverse table has been extended from five to 
six tabulated decimal digits of kK for 


V = 0.01(0.01) 0.99 


(14) 
ec = 0(0,01) 1.00 
and the table has also been supplemented with an extension of the 
variable | from 0.99 to 0,999999 such that 
1 = .9906(.0005) .9960( .0001).9999(.00001).99999(.000001) .999999 
(14a) 


e = 0, ¢ = 0.10(0.05) 100 


This table is given in Appendix D with a corresponding inverse table 
for the circular coverage function, P(Xk,d) in Appendix E, This P(k, d) 
table is given by a tabulation of R as a function of P andd. Poth 
tables are correct to within one unit in the last digit given. The 
ranges of / and d are as follows: 


P = 0.01(0,01) 0.99 
(15) 


d = 0(0.1)5(0. 2) 10( 2) 20( 5) 120 
and 


P = .9900(. 0005). 9990(, 0001) .9999( . 00001) . 99999( . 000001). 999999 
d = 0,.05,.10,. 25, .50,.75, 1, 1.5, 2, 3, 4, 5, 6, 8, 10, 20, 30, 50, 80, 120 


15a) 
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Recently, L. Harter [6] published a table of I'(k,c) versus K 
and c and an associated small inverse table of K as a function of | and 
ec. His approach however is limited since the quadrature method he 
uses breaks down or requires a prohibitively large number of inte- 
gration intervals for values of ¢ near zero. Furthermore it is not 
efficient becausz it is not known a priori how many integration inter - 
vals will be required for a given input. Neither of these deficiencies 
is present in the program to he described. 


The equations described for Il(k,c) had been worked out 
a few years ago by one of the authors of [11] and independently by 
others, [2], [3]. The success of the method for computing I'(k,c) by 
these equations warranted considerations of extending the technique to 
the P(A,d) function. This was easily possible because of the existence 
of a functional rclationship between P(k,d) and I(A,c) established in 
equation (18) of the next section. This relationship is derived 
by H. E. Fettis in reference [3]; however its use was 
reported by him for a somewhat different purpose than to develop an 
improved technique for computing P?(h,d). The merit of the extended 
analysis included herein is the greatly increased speed and accuracy 
attainable in the computation of the circular coverage function, P(A, d). 
relative {o other existing methods, by the use of Fettis' relation, The 
capable speed and accuracy of the computing program for P(h,d), as 
wellas I'(K,c), is manifested by the two inverse tables mentioned 
above. Mach table required the computation of approximately 45, 000 
P(k,d), or V(K,e), functions to an accuracy of seven or more decimal 
digits. The NORC computing time for each table was less ihan cight 
minutes. 


IV. RELATIONSHIP BETWEEN P(A, d) ANDI'(K, ¢) 


The relationship between ? and | can be derived by utilizing 
two of Fettis' preliminary results which are given in terms of g = 1-P 
in equations (I - 35) and (I-44) in [3], They are given in terms of P as: 


; ja -d| (4)Rx~d J 
(Ad) —P(dh)yaa VF kdl, 16) 
P(Kd)—P(dh) = 4 ( Nae (22 ‘ 

_ Red? 
P(K,d)+P(dh)=l—e 2 1, (Rd) (17) 
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Equation (17) is easily derived; the origin of equation (16), which is 
given by Fettis in a slightly different notation, is unknown to the 
authors. 


It follows, by adding equations (16) and (17), that * 


K? +d? 


ahd f. |A-d|\) (4) A >d 
l-e 2 (kd) 4b {IK dj, (18) 


P(h,d) = 


nile 


Rad (-)R<d 


Thus a value of the relationship, as given by equation (18), is that if an 
efficient computing program exists for computing the !'(k,c) function, 
then an efficient computing program can also be realized for computing 
the circular coverage function. The basic features of a computing pro- 
gram will be presented in Section VI whereby | (k,¢) is computed © 
directly and ?(k,d) with the aid of equation (18), 


Y. RECURRENCE RELATIONS FOR 1'(4,c) AND al’/dk 


The function dl /dk is necessary for computing the inverse 
function A(I,«) by a Newton-Raphson procedure (See Appendix C), The 
function 0! /dA canalso be conveniently computed simultaneously with 
1; therefore relations for both functions, | and d/'/dK, are developed 
in this section. 


A few simple linear transformations and the substitution indi- 
cated by equation (11) reduce the integral of equation (4) to 


ARR /4 
V(A,e¢) = = exp [-27 "| 1 (20) de (19) 
0 
Then 
(x) ON ee ge BK 2p CaKty2) (20) 
BRJ OK 14 c¢* 2 


Guenther recently (see equation (2) in [16] ) derived an equation 
for P(A, d) in terms of I(x) and the incomplete gamma function, which 
can be shown to be equivalent to equation (18) of the present report. 
However, he did not exploit his relationship from the point of view of 
developing an efficient program for a high speed digital computer. 
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Furthermore, the following relation 


: 2 [ : ] T Ce) d : ! 1 
—_— exp f- —u SS = 
Ae A e Ac KK RB / A y? ery (21) 
oe 


(see page 76, Example 14, of [5]) is used, with a linear transfor - 
mation on w, to write equation (19) in the following form 


i om ed B tw 
V(RK,c) = 1 - ai [- x] is) du (22) 


The properties of |,,(*) are classical and are exemplified by the series 


expansion definition: 
oo 1 2 x 2n 
Kio (3) (3) (23) 


n=O 


which is uniformly convergent on any finite interval of x. A semi- 
convergent or asymptotic series expansion for [,(x) is given by: 


NV 
exp (x) [(2n) ']? ox yn (24) 


I (x)~ eaLeDES eine 
o V 20x Par BAN (Cn lys 


i,e., the right hand side of equation (24) represents |,(x) only for 
sufficiently large values of x and finite ¥, The extremum values of x 
and V are determined whenever the precision desired in I,(*) is speci- 
fied, (See Appendix (A) for further comment), 


The integration of equations (19) and (22) is accomplished by 
substituting the series given by equations (23) and (24) respectively, 
with proper substitutions forx, in terms of «, and subsequently inte- 
grating each resulting series term by term onuw. Thus 


pee gang thle 5 co 
V(A,e) = re > tall exp (23") ww? diye = De Te (25) 
° 


n=O 


N ve N 
oe I [( Qn) 1]? 1 
V(K,e)~1- = ——— | exp ( x") we (2041) dw =] = > df 
Zoya go, 2a) mA a 
4 


nF AK? n=o (26) 
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The term by term integration of the series which results in 
equation (25) is justified by application of the Weierstrass ''M" test, 
since the series before integration is bounded by 1, (4K*/2) for all 


finite values of 4K?/2, i.e., , 

~ — 1\? LAN PAREN" AR? 

be f, Ces > exp (~2Bw/A) ea wer ¢< > (=) e =I, =o (27) 
n=-o n=0 n=O 


for all w on [0,AK?/4] 


The term by term integration resulting in equation (26) is also justi- 
fied for all values of AK? for which the relation (22) is valid, because 
of the existence of the integral given in equation (22), (See page 17 of 


[1]). 


The explicit computation of each term, after the first, in the 
series of equation (25) is achieved rapidly and accurately through the 
use of recurrence relations, Let 


Iv 


“‘ 


AK? /4 ‘ 
= | exp ar w) ww?” dw n>o (28) 
o 


A repeated integration by parts leads to 


4 A MK? nA (AKA! ee 
tn B Varteee EGET) eM a2 ca 


The terms 7,,, that appear in equation (25), can now be identified as 


a 2 1\2 . 
ben Ney) Yen (30) 


If the function |,(1K?/2) is replaced by its series representation as 
given by equation (23), the individual terms of the resulting series in 
equation (20) are given by 


: AK?\? | 1 f 1\2/AK7\?" 3 
Sep = e ) Sapa = ry (=) (Z) exp (-BK /2) n>] (31) 
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2c 2 
5, = manga exp (—BK*/2) (32) 
lee 
where 
e 2 a 2 
x , e exp (-BK?/2) (33) 
»> eR l+c? ° 2 


Substituting the expression for ¢,, from equation (29) into equation (30) 
and using equation (31), it can be shown that 


Qn ~ 1 1~— ¢?\? Aa 
Yan oF. ( \ ( ) Cages (1 “ayn n> (34) 
MN ] "2 ae, } 2n 2 pie 2n 
‘ / \ 


GA 
oN 


The initial term for» = 0 follows directly from equations (28) and (30) 


9 » : 
Sofi ieep tank? 7 21S ; : z~— 0 (35) 


1+ c? +e 


To = 


oO 


One may observe that if the sum given by equation (33) is accumulated 


as one computes | then 
ov we 
oes > en (36) 


n=O 


In the case « - 1, which was mentioned in Section II, 4 = 0 and 
k= 1 by equations (6), Then 4, = exp(-A?/2) and 1, = 1 - exp(-A?/2), 
by equations (32) and (35), It is seen from equation (31) that all the 
s’s except §, vanish (since .! - 0) and thenfrom equation (34) that all the 
7’s except 7, vanish (since | -c*? ~ 0), Equations (25) and (36) then give 
M(K, 1) = 1, = 1 exp (-K?/2) - OV(K,1)/0K = K exp (-K? /2) 
respectively, which are the correct values for this case (c = 1), 
‘Thus no modification of the computing program is required for this 
case, 


The explicit evaluation of each term of the series in equation 
(26) can also be accomplished efficiently by suitable recurrence 
relations. Let 


Longs 7 |» (-0 / 2A) wo A) die (n > 0) (37) 
AK? 


10 
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An integration by parts yields 


2 “5 —(22=1) Refs 1 ; 
Loney = Onl (ARS) ? ¢, anit By Fant) (n > 1) (38) 
By incorporating equation (38) with equation (26), it follows that 
5 1}2 
‘ __! [ (an) 1] one 
ant 2hefa DANE 153 ane n 20) (39) 
where 
ft o 
L,= | exp (= /24) wT! ? dw (40) 
AK? 
The substitution 
ya few, 2 (41) 


transforms the integral in equation (40) to 


Ly bye | esp (=y? dds (42) 
EN fe 


If one considers the equation 


ay ee ys fin? ; 
trae) = AC eet BA? 72 u() (43) 
C Cc 


AK 
and replaces the function [, “*) by its asymptotic expansion as given 
by equation (24), then the individual terms of the resulting serics for 


silicone 
(l-c Ure are given by 


4 
‘ | c\? an — | . 
eee - ( _ =) @) ( n ) Makes (n 2 1) (44) 


where 


11 
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and 


ov 
is (46) 


Furthermore, it is advantageous from a computational standpoint to 
introduce another function, Y,,_,, which is related to X,,,, by the 
following: 


Te? [( 2)! — 
x =(% -1)Y = pote Uae) HT" (««*) 2 exp (—K?/2) (n 2 1) 
Vi 


anti an-1 Vv 2 Qtr nt) 
(47) 


1 c¥ fx -1 
er a i =) (z) (7) Renn (n > 1) (48) 


Replacing .,,,, and /,,_, in equation (38) by their expressions in 
terms of the corresponding ¥,,,, and ¥,,_, from equation (39), sim~- 
plifying, and using equation (47) leads to the following recurrence 
relation: 

eee ees y <A 

an+1 _ a 7 On “anal ne (49) 


—C hake 


such that 


If equations (39) and (42) are combined, an initial term (n = 1) of 
equation (49) is given by 


oo * 
a : exp (-y¥?) dy; (50) 
‘Vine VF af 
K/f2 
* The initial term 4, as given in equation (50) can be written as 
My = jee 1 ~ Frf(K/¥2) | 


i lye 


The computation of the £rf function and its derivative, given by 


2 2 
Erf’(x) =s—e*, 


Va 


on NORC, is done by subroutines designed by Dr. A. V. Hershey which 
restrict the error to the thirteenth digit. A faster method for computing 
Erf(x) and Erf’ (x) ata sacrifice in available internal machine storage 


12 
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the other initial term of equation (49) is obtained by'setting n equal to 
one in equation (48) so that 


: 1 c \2 ] 
= mate (z) (7) Xx, (51) 


The necessary recurrence relations that will be required for 
computing P(k,d) and aP /ak from equation (18) and (13) respectively 
can be derived from those already given for | (A,c) by appropriate 
interpretation of the variables. This will be siiown explicitly in the 
next section wherein the order in which the relations are used will be 
described for P(A, d), AP /AR as wellas I'(K,¢), AV /OK, 


Vl. COMPUTING PROGRAM FOR !'(k,c¢) AND P(R, d) 
A) Program A 
1, Operation: 
Program A is used to compute the probability function 
V(A,c) from given values of A > 0 andO<c <1, Simultaneously with 
the computation of !’ the function @!'/dK is also computed, although it 
is not available as a direct output quantity of the routine, 
2, Formulation: 
I'(k,e) is computed from equation (25) and 0J'/dK from 


equation (36) if 
AK? < I 


locations is given in Appendix B, The method of Appendix Bis used in 


conjunction with the IBM 7090 program for computing I and P, If 
e= 0, K > 0, this equation for ',gives "= 1 &rf(K /¥2) and itis 
seen from equation (48) that ), and all subsequent }’s vanish, and 
hence from equation (49) that all the 1s after !, also vanish, Then 
equation (26) gives I'(K,0) = 1- V, = Erf(K/¥2). This justifies the 
statement in Section IJ that I'(K,0), if defined as Erf(K /J/2), is com- 
puted by the regular procedure and need not be regarded as a special 
case, 
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if AK? >M ” 


then (K,c) and dV /dK are computed from equations (26) and (46) 
respectively. The choice of 1! is discussed in Appendix (A). 


If AK? <M 


the evaluation of equations (25) and (36) is started with the following 
initial values: 
4 
n= 0 > = 0 d= 0 (52) 


2 2c 
exp (-BK?/2) T, = [ - 5 | (53) 


S.=s— 
oO 
L + c? 


The steps in the iteration are given by 


Qn + 2—» Qn (54) 


2 
AK? ; 
( An ) Sane en (55) 


rs a ==) : 2 eee (56) 
eer a an—2 ~ ones San lan 
rE (51 
ra 7 
> Sin >), (58) 


The iteration is recycled until 


S.,$6= lo 
(59) 

Ton S 6 = 1077 
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where jis a positive number. The final results are’ given by 


K(K,e)~T, +>, (60) 
/ 
av /8K ~ BK [s, > | (61) 
If AK? > It (611) 


the evaluation of equations (26) and (46) is started with the following 
initial values: 


n=0 >> = 0 <0 62) 


: L~c? 2 °° 
Ay = a exp (-K? /2), Moos poeta ee exp (—¥?) dy (63) 
2 va ' 1 -c? vn pt) ; 
K/¥f2 


The steps in the iteration are given by the following recurrence rela- 
tions: 


Qn + 2—» Qn (64) 
| c\2 fn -\ ; 
Tap? kK = in Me tee akin (65) 
K . ? Ons 2] i 
ae eed ee ea (66) 
ae ig 


(2 — Way am Na (67) 


y+ Wonat7de (68) 


éf 


/ 
> a eran Fay (69) 


The iteration is recycled until 


AG eees 
(70) 
V <e = 1077 


2H = 


15 
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The probability function and its derivative are finally given by 


V(K,e)w(1-¥,) = (71) 


av /an~ ( 


) [x + (72) 


1— c? 


B) Program B, 


1. Operation: 


Program B is used to compute the probability function 
P(k,d),and dP/dk from given input values of A > 0, d 20. 


2. Formulation: 
PCR, d), ok are computed from equation (18) and (13) 
OR 


respectively. The recurrence relations for the quantities | (xd), 


\ (Via i =t) which appear are easily obtained by the following 


R+d 
substitutions: 
|A — d|-»K (73) 
Rd 
tay ara (74) 
KR+d 
Then if 
QRd—» AK? <M (741) 


the evaluation of equations (25), (33), (18), (13) is started with the fol- 


lowing initial values 


n=0 Y= 0 Ye 0 (75) 
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im rR? y d?\ * R? — d?| = 
amas ( 2 ) mer {#=9) b- | ae 


The steps in the iteration procedure are given by 


Qn + 2» 2Qn (77) 
kd\? . oa 
(=) Sanna an n>1 (78) 
Rd fx 1\ |R? ~ d?] ; An : : 
aa es F 
Re 4 d? mn "ees hk? 4 d? R? + d? an * “an (79) 


>: ono > (80) 
en Son x (81) 


(82) 


* If the variable $,, of equation (31) is expressed in terms of & and 


t 


d by means of the equivalences indicated in relations (73) and (74), the 
F [ak? — d? | Re + d? ] 2 kd 2n 
result is Son =o exp {-— — —j} - In evaluating 
Re 4d? 2 n} 2 


P(R,d), however, it is convenient to omit the factor |A? ~- d?|/(k? + d?), 


which is not present inthe expression for o6P/dK , equation (13), The 


sy k2 = d? _ 
resulting variable is denoted by 5,,. Thus S, = pie so 
R? 4 d? 
7,, has the same significance here as previously (Section V), and 


2n 


hence no new variable hy is introduced. 


17 
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The circular coverage function, P(R,d), and its derivative with respect 


tok, i.e., o@P/dR , are given finally by 
1 ~ (+)R>d 
Rego: |e [25 x 
ener Al (5, )e(r >) crea &? 
oP oa ‘ 
pare ; (84) 
oe” * (+5) 


~ R—d 
Since > Ti, =F (in - dj, E 7) < 1 (see equations (18) and (83)), 


n=O 
and since each 7,, is positive, as can be seen from equation (25) if 


interpreted in terms of X and ¢ rather than K and c, it follows that each 
individual term 7,, is less than unity. It can also be shown from rela- 
tion (79) that the 7’s form a strictly decreasing monotonic sequence, 
each (except 7, ) being formed from its immediate predecessor by 


multiplication by a positive constant (supposing d > 9) less than unity, 
followed by subtraction of a positive quantity. If d=0, the P(k,d) 
function is evaluated analytically as pointed out in Section IJ. The 


largest T,, is 7, 


= rc k? a qd? ? L ) a7\2 > y 
>> Son = exp [- 3 1, (Ad) = exp | — 3 (kK -w) exp (~Ad)1, (Ad) 


n=G 


as given by equations (76). Since 


(also by comparing equations (18) and (83)), it is easily shown that 


a. 1 
Ds San < &XP - ae as 1 again,assuming that d>0. For the 


nh=o 
function f(x) = exp(—*)I,(*) by standard properties of Bessel functions 
has a derivative — exp(—)[I,(x) - I,(+)] which is negative for all 


values of x, and hence,for positive values of x (or kd), this function 
has a smaller value than its value at x=0 , whichis unity. Fach 


18 
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bs rea d\ (1 \* (ra\ | 
oo ei a a ny Vo is positive. If kd <2, itis clear 


from the recurrence relation (78) that the largest of the S’s is oo 


which is given by the first of equations (76). If Ad>2, the S’s 


increase to a maximum term and then decrease. It can be shown that 


1 Rd 
the largest term does not exceed exp @ A (h -a)/ (. ze \ ines 


R F : td 
[42] is the greatest integer not exceeding ss 


yu Okd > i (841) 


the evaluation of equations (26), (18), (13) is started with the following 


initial values: 
/ 
n=0 > = Q = 0 (85) 


y ya wis ] I K+d\ 2 7 
'< iy) Pamrre ee Pa 2 su 
x, v2kd exp ( mt ar) ; W, ” ard (Aes ‘- exp(~y”) dy (86) 
|K=d| //2 
The steps in the iteration are given by the following recurrence rela- 
tions: 
Mn + 2—» Qn (87) 


x here does not correspond exactly to ¥, of equation (45) under 
the equivalences indicated in relations (73) and (74), and hence the 
change in notation. The relationship is ¥, = (k 4 d)¥ 


¥ 


and generally 


=(R+d)X Similarly a new variable ¥,,_, =(K+d)¥,,_ 


en+it 2n+ ft i 


1s used in relations (88), (89), (90). See an earlier footnote regarding 


of the P(R,d) program. I 


e,4, nas the same signifi- 


the variable a 


cance here as previously (Section V), and hence no new variable 


Vena, is introduced. 
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- 1 {m—1\ — ph 
‘Ikd ““2n—t | heme {88) 
12 ee {? a k = r.2 Oy Z. l 
ES tcl me ee Gee We sees (89) 
AR ant Ap 9 an— andi 
tRd td Nn 
es 0 ee ree (90) 
ys Und, (91) 
é 7 


The iteration is recycled until 


Neh Rie (em 
F (93) 
Yon 21 SES lov! 
° The circular coverage function, piyx.d), and its derivative with respect 


to h, dPoah , are finally given by 


MCL ad ste [he] So eee ey ' [ns oe 
sees ;| vai ee I | (-)R <d 
mo, MSS ‘ 


In the case Kk -d, which was mentioned in Section II, equation 


(18) gives Pth,dy - (172) [1 = exp(-k7)1,(h7)) . It is desired to show 


that the program gives this result without modifications. The quantity 
Rd, or 2h?, could satisfy either of inequalities (74 I) and (84 1). In 
the former event, it is necessary merely to observe that all the 7’s 
vanish, by equations (76) and (79), The correct values of P(A,d) and 

oP dak are then given by equations (83) and (84). If 2d = 2k? >, 
satisfying inequality (84 1), , - 1 by equations (86), since 


(2/¥m) | exp (-y7) dy= Erf(o) = 1,and equation (89) shows that all other I's 
vanish.” P(k,djand dP /ak are then given correctly by equations (94) 


and (95). 
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Vil. DISCUSSION OF THE TABLES 


The tables are identified as Table I and Table II. Table I con- 
tains a tabulation of K as a function of V and c. There are two ranges 
of and ¢ in the table. The first part of the table has the following 
ranges for lV’ and c: 


—< 


’ = 0,01(0.01) 0.99 


(96) 
= 0( 0,01) 1.00 


a 


econd part of Table I has K tabulated for the following values of 
ce 


Ve = 0,99( ,0005 ) 0.999( .0001 ) 0,9999( .00001 ) 0.99999 (0.000001 ) 0.999999 (97) 
cz 0, c= 0.10(.05) 1.00 


The values of A were determined by using a Newton-Raphson 
procedure as described in equation (149) in Appendix C. The value of 
|’ for which K was desired, ?, was held fixed and a sequence of values 
of K were determined by letting ¢ range through the increasing values 
of « indicated above. The value of P was then incremented and the e 
values were spanned again. Thus the order of the computation con- 
sisted of computing values of K across each horizontal line of the table 
and then proceeding to the next horizontal line in the direction of in- 
creasing P. 


‘The Newton-Raphson procedure requires a starting value of A. 
Each starting value of K was determined by choosing the final value of 
A computed for the previous value of c¢, except for ¢ equal to zero. 
The starting value of K for « equal to zero was chosen as the value of 
kK computed for the last value of P with ¢ equal to zero. 


The value of K was accepied as correct to within one unit in the 
sixth decimal digit for any P and ¢ in the range specified by the 
equalities (96), (97) whenever 

V(K,,¢) - P| <1 x 107? 
| Vik, | x (98) 


(Ke Reo hes 2107? 
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‘The value of j, that appears in equations (59), (70), and '! were chosen 
such that 
j = 8, df = 30 (99) 


for P and ¢ satisfying equalities (96), and 
j= 12 i= 40 (100) 


for P and c satisfying equalities (97), The values of I! were chosen 
conservatively with respect to the analysis given in Appendix (A), An 
average of four iterations were required to determine A for the range 
given by equalilies (96) and an average of five iterations for the range 
given by equalities (97), The average time of computation per case 
is ten milliseconds for NORC; this time could be reduced to about six 
milliseconds which is the average computation time on the IBM 7090 if 
Erf(a) and Lr{‘() were computed as on the 7090 and as described in 
Appendix (#3). 


‘Table Il contains a tabulation of & as a function of ? and. There 
are two ranges of ? and d. The first part of the table has the following 
ranges for ? and d: 

P = 0,01(,01) 0,99 
d = 0(.10)5(.20) 10(2) 20(5) 120 cron) 
‘The second portion of the table has & tabulated for the following values 


of P andd: 


P = 0,99( 0005 10.999 (0001 ) 0.9999 (00001 ) 0.99999 ( 000001) 0.999999 (102) 
d = 0, 0.05, 0.10, 0.25, 0.50, 0.75, 1.0, 1.5, 2, 3, 


4,5, 6, 8, 10, 20, 30, 50, 80. 120. 


‘The values of k were determined by using a Newton-Raphson pro- 
cedure as described by equation (150) in Appendix (C). The value of 
P(R,d) for which A was desired, P, was held fixed and a sequence of 
values of &k determined by letting ¢ range through the increasing values 
indicated above. The value of P was then incremented and the /¢ values 
were spanned again, The order of the computation consisted of com- 
puting values of & across each horizontal line of the table and then 


proceeding to the next horizontal line in the direction of increasing P. 
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The Newton-Raphson procedure requires a starting value of R. 
This starting value was chosen as the final value of 8, which was tabu- 
lated, for the previous value of d, except for ¢d equal to zero, The 
starting value cf Kk, for d equal to zero, was chosen equal to the value 
of x tabulated for the last input value of P when d equalled zero. 


The value of & was accepted as correct to within one unit in the 
last digit position given for any P and d in the range specified by 
equalities (101), (102), whenever 


1P(K,,d) ~ P| <1 * 1077 (103) 
hy ate he S10? (R, <4) (104) 
ay ieee mine lee al ceo a (R, > 4) (105) 


n 


The value of j, that appears in equations (82), (93), and !! were chosen 
such that 


i=8 Ve = 30 (108) 


for P and d satisfying equalities (101), and 


oe Wo = 40 (107) 


for ? andd satisfying equalities (102), The values of ¥ were chosen 
conservatively with respect to the analysis in Appendix (A), The num- 
ber of iterations and computing time are identical to those already 
stated above for the |’ function. 


After the radii xk and & had been computed, by the Newton - 
Raphson procedure, satisfying inequalities (98) (for x) and (103) - (108) 
(for Rk), they were checked by o method described and flow-charted in 
Appendix (C). 
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rn A St 


ESTIMATION OF THE FACTOR M 


This section is devoted to a method for estimating the lower 
bound of ¥ that could be used in the inequalities (611), (741). A sub- 
sequent conservative estimate is obtained for the maximum number of 
terms required for the Taylor or asymptotic representation of F(K,c) . 
It will be indicated that the value of M chosen for computing }(K,c) 
will also be satisfactory for P(k,d) at the same ievel of accuracy. 


The estimate for !! is based on two requirements which are as 
follows: 


a) at least one term of the asymptotic series of equation (26) 
must attain a value less than a pre-assigned «>0, for all 
K ande 


whenever 
AK? > (108) 


b) the Taylor series and asymptotic series should be trun- 
cated at approximately the same number of terms in the 
neighborhood of 

AK? = ¥ (109) 


Condition (b) can be relaxed somewhat since it is not a critical re- 
quirement. It is imposed primarily to maintain the computing pro- 
gram at more or less optimum efficiency. This is based on the 
assumption that the computation time per cycle as well as the number 
of exponentials, rf functions and square roots to be computed for 
each series is equal. 


Condition (a) however must be satisfied because the final results, 


for AK? >, depend directly on that requirement. Furthermore, even 
though one or more terms of the asymptotic series, all of whose terms 
are positive, do become less than epsilon (-), the stipulation that 
V(K,c) will be computed with a pre-assigned accuracy, ¢ , is not as- 
sured a priori, The reason for this is that no rigorous error bound 
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exists, in general, whenever all terms cf an asymptotic series have 
like sign. Thus in the final analysis it was necessary to check the 
results, derived below, by numerical experiment. The experiment 
is described in this appendix following the analysis. 


The procedure for estimating ¥ is an iterative one, i.e., various 
values of AK? are tested until condition (a) is satisfied; simultaneously 
the asymptotic term of minimum value, say the WV, th term, is deter- 
mined. Upon the determination of \,, the value of AK? is substituted 


into the corresponding Taylor series, equation (25), and the first term 
of that series, which becomes less than the epsilon, say the 


Nyt term, is determined. If 
[Ng -Npl >] (110) 


the procedure is iterated by changing the value of AK* in such a way 
that condition (a) continues to be satisfied and the magnitude of 


IN, -Nr| is reduced. 

In order to be specific, let «= 1x 10% and assume six correct 
decimal digits are required for the final result, (K,c). Then N, is 
determined in the following way. The general nth term of the asymp- 
totic series, equation (26), is given by 

1 T(2n)1]? 
Manas = EN 2" (n1)8 [on-esane (2444) dw (n>0) (111) 


Ak? 
Then 


l On) 1 |? es anti 
M < [Gea ent-a*/2f ah 2 ) tw 


Ut 


an+t QAeVa 24" (n1)3 


2 2 
d ieee (x2 /2) AK bee (lean <N) (112) 
= Qhefn 24"(n1)? (2n-—1) 
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The use of Stirling's approximation for the factorial in the preceding 
inequality leads to 


M 5 ese eke ao (113) 
an+1 oa aes Gea Ge (n 2 1) 


In order to satisfy condition (a), for «=1%x 10-*, it is apparent that 


AK? must certainly be greater than ten. Moreover if six correct digits 
are desired in the final result, then since for a fixed IV (K,c) the radius 
K is amaximum fore = 1, it follows directly by evaluating K from the 
equation 

V(K,1) = {1 - exp(-k*/2)] = 0,9990905 (114) 


that the range of A, under these conditions,is limited to 
O<K < 5.26 (115) 


For values of AA? that satisfy the inequality 


AK? > 16 (116) 
the range of c’ is limited by the inequality 
O<c? < 0.47 (117) 
Using the above inequalities, it follows that 
} : < 1.2 
ae : 118 
Qn V1 -c? ( ) 
Therefore 
Mees < 12 aes | z eat) (119) 
Vn in=1/2){ | tke 


By considering f(») as a differentiable function of 2, the approximate 
minimum value of I/,,,, is obtained by determining the root of the 


equation 
‘ f (ag)/flg) ~ 0 (120) 
The implicit relation for 2, as a function of Ak? is 
1 2 
AK* = 1, exp | - — - ——— (121) 
| a, 2,1 
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The substitution of this relation into the equation for f(n) gives 


' Oa ee 1 Oh, 
WM spa SE ie on - exp barre (122) 
Dropping terms of 0(1/n,) gives 
2 
AK"~ ng (123) 
= 3/2 3/2 / 
Mpg pee boty e 0 O12 e SAK2 ) os goin (124) 


Thus for 
AK? = 16.25 


N4 ~ 17.0 (125) 
Wana , (min ) <107° 
The procedure continues in order to determine whether 
Ve~ 17.0 
such that condition (b) is satisfied. 


The A term of the Taylor series is given by equation (25) as 
AK? / 4 


_ ae Ae Ue 
toy ot rr @ [os (=2Bw JA) we?” duc (rez 0) (126) 


The use of Stirling's approximation for the factorial and the first mean 
value theorem of integral calculus implies 


2 1 e an Ak? / 4 
ye 7 20 5 
Lon < le Dan (:) 02" de (127) 
oO 


de tps | Re Ke AT 
ee e 9) 2n + -( n ) ey An ) Cee bee 
(128) 
K? l dk2e “| 7 1 [Ake m 
the ( An ) ae ( ) bee Ses t 
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For 


AK? ~ 16.25 (129) 
Noel ea (130) 
heat 2.23 (ee). ee ae 
“Taz? \ 68 
whereas 


T,~ 6.0 * 1078 
Thus conditions (a) and (b) are satisfied te io?) by 
Vl 16,25 
(132) 
Vr Ny Iv 
In order to carry oul a similar analysis for the Pufl./: fune- 
tion, it is necessary to consider the second term of eguation (18). Tf 


the asymptotic series is used to represent 1, then the a” tern, of 


oO? 
the series which represents the second tevin in equation (18) would be 
| fy exp [- ( R-d za] (2 ty? as 
—-- ee ee - oe a Ta 
Wd 2M J2akd 


In the case of the ‘Vaylor series, the no term would be piven by 


; : | Aid oH Q. 
noe exp - (n? od?) | i) ( ) (134) 


vy 
After introducing Stirling's approximation and replacing the exponential 
term by unity, the quantities Wand . are aporoximated from the 
relations (133) and (134) as 


te Rd ett — l nM se ie ine 
eae a ame ear oy ce ae 
*N tan \ HW Wd 2 Qnkd Ye \2Rde 


‘Tne substitution of unity for the exponential terms leads to conservative 


(133) 


a: 0 


estimates for If and V¥. The value of Wr - 7") was 
ORs AP IGS, (136) 


and subsequently, the value of \ was found to be 
Nee hrs 
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ea A AE SIT PAS A A a SS A ya r eNAE 


The analysis above is crude and is given only to indicate the ap- 
proximate values of i and W that one might choose if all other factors 
were equal,i.e., if the asymptotic and Taylor series required computa- 
tion of the same number of exponentials, square roots and £rf functions 
and if the number of operations in the basic computing loop to calculate 
an arbitrary term of each series were the same. Actually, this is not 
so, even though the loop is about the same length in both cases, because 
use of the asymptotic series requires the computation of an Frf function, 
exponential function and square root whereas the Taylor series requires 
the calculation of two exponential functions for the IBM 7090 program 
and only one for the NORC program. Thus in order to maintain effici- 
ency and also insure sufficient accuracy in the use of the asymptotic 
series, an " factor equal to thirty was chosen for the NORC program 
and an "4 "'' of twenty for the IBM 7090 program. The maximum num- 
ber of ‘Taylor series terms used in the NORC program such that the 


series were terminated when any term became less than 10-* in value 
was twenty and about sixteen for the IBM 7090 program. 

The computed results were verified as correct, to within the 
specified epsilon, by the following experiment: 

A range of values of 4 and kK were chosen which spanned the do- 
main of v(or?) in the neighborhood of AK? (or 2Rd ) = 18~30. Then 
each individual case was run on the computer using first the Taylor 
series and then the asymptotic series to compute the value of (or P) 
for the same value of K (ork), « (or d). In all cases the values of 
Vior P ), by this "overlapping" technique, agreed to within the accuracy 
of the specified epsilon. The minimum « hada value of lo’. It was 
concluded therefore on the basis of thesc tcsts that the series given 
herein would yield results of prespecified accuracy on the basis of the 
criteria given above in spite of the lack of rigorous practical error 
bounds on either series. 
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APPENDIX B 
COMPUTATION OF THE ERROR FUNCTION (#7/(x) ) AND ITS DERIVATIVE 
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COMPUTATION OF THE ERROR FUNCTION (érf(x)) AND ITS DERIVATIVE 


x 

2 

The &rf function, Erf(x) zal et? dt, (137) 
0 


2 2 
and its derivative, Erf (x) Bey e~ , are computed rapidly and with 

7 
sufficient accuracy in the present IBM 7090 program by using an iter- 
ative interpolation procedure as described in [7]. See also Sec. VI of 


{12], from which the following account is taken, 


The functions, #rf(x) and Erf’(x) have the following Taylor 
expansions: 


K K 
( Ax)" 
Bria) ~ YW) 4) = De ty (138) 
neo Me n=0 
K K 
a, (Axy 
rf! (2) SY BO) ag) Yoon (139) 
n=o n=O 
t 
where Ax =x -%, = py such that |p| s= . 


The quantity «+, is representative of the values of « for which 
‘Erf(*) and Erf’(x) are stored at intervals of length y such that 


1 
lx ~ #91 S57: (140) 


The functions ”'") (x) for a z 2 are equal to the classical Hermite 
polynomials, //,_, (x), multiplied by Erf'(x), [7]. In defining the 
Hermite polynomials, some works on special functions use different 
conventions regarding algebraic signs. The present authors follow the 
conventions of page v in the Introduction of [7], wherein 

Hy(%)=-2, I(x) = 427 ~ 2, (x)= —8x° + 12x, .... The recur- 
rence relations which generate the individual terms /,, L, are given as 
follows: 


NWL REPORT NO. 1768 


Let 
In = Tht (141) 
Loe FW), Ge 8,84. oe RG (142) 
where 
Ax 
In am Iya, Jy oe (143) 
HO?) (x)2—2 [we-? (x) + (n-2) 4-2) | (144) 
Thus 
> At Ax 
ce tare E Lest (5) (n~2) ea (145) 
Ax 
L, =- 2 mt [soln + Nas | (146) 
The starting values are given by 
I = Wag) = Erflxy) 1, = Ard!’ (x,) = AxErf (x,) (147) 
/ .,f a / 
Lait (a) = Erf (x,) LL, = Axli (xg) = = 2x, Axil(x,) (148) 


If y is sufficiently small, very few terms of the series given by 
equations (138), (139) are required to evaluate the Erf function and its 
derivative to a specified accuracy. The integer « is therefore a function 
of epsilon and gamma. It was determined for different e and y by com- 
puting the maximum value of u'™) (x) on the domain 0 S$ %_ § 3.85, (see 
equation (1) in [7] for error term involved). This was sufficient to 
bound the absolute truncation error. A table is given below for the 
absolute truncation error as a function «x andy. 
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je] y= 0. 10 y = 0.05 


3,02(~4) |7.56(—5) 
5,88(-6) | 7.35(—7) 
7.17(-38) | 4.48(-9) 
3.44(-—11) | 3.53(~13) 
1,25(-11)] 1.96(-13)|8.01(— 16) 


Absolute Truncation Error for Various Values of «x andy. 


Values in the Parentheses Indicate the Power of 10 by which the 
Tabulated Values Should be Multiplied. 
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APPENDIX C 


PROCEDURE FOR COMPUTING AND CHECKING TABLES 
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PROCEDURE FOR COMPUTING AND CHECKING TABLES 


In table{, K is determined asa function of V andec. This is 
done by the Newton - Raphson procedure. Thus the n iterate of 
K, K,, for a fixed value of Y, equal to P, is 


flee os, eee, n 
n not av(K,_ |. c)/ 0K 
Similarly, in table Il, k is determined as a function of » andd. 
This is also done by the Newton - Raphson procedure. _Thus the nth 
iterate of R, X,, for a fixed value of P(k,d)', equal to P, is 


1 (149) 


Wy 


PRS 225 2G P : 
Pp = Rf a sitet n> 
Bn Kaas OP (R d)/9R 7 ee 

All values in both inverse tables were checked, and ina few 
borderline cases corrected by one unit in the last figure retained, bya 
checking program basec solely on the direct computation of probabil- 
ities. The details of this program are shown in the accompanying flow 
chart, Although the letter 4 is used in the flow chart to denote the 
radius of the circle corresponding to the specified probability P, this 
procedure was used for checking both the V’(k,c) and the P(R,d) inverse 
tables, The last approximation for the radius (Kors) determined by 
the Newton-Raphson procedure is denoted in the flow chart as R,- This 
value &, is then rounded toa predetermined number of figures, seven 
significant figures if &, - 1 and six decimal places if k, <1, and the 
rounded number is called A, and is the value of the radius given in the 
inverse table if subsequently verified by the checking program. The 
essence of the checking method can now be described in one sentence: 
w= P(R,~h/2,d)-P and v=?P(K,+h/2,d)~P are computed, and, if 
ug0 and »20, k, is thereby determined as the rounded value of the 
radius, correct to the last figure retained, The symbol 4 here repre- 
sents the value of one unit in the last place retained; for example, if 
10 <¢R,< 100 so that Rk, is carried to five decimal places, then 4 = 107°, 
If the conditions u <0 and » 2 0 are satisfied, then Rk, is bounded be- 
tween &, -h/2, which is too small to give the specified probability P 
(unless u = 0), and XR, +4/2, which is too large (unless » = 0), assuming 


nowt? 
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that the direct program computes probabilities correctly, But ifu or v 
is excessively small numerically, say less than 5x 1l07'', then there 

is doubt as to the correctness of the last figure evenifu<0o, andv>4@, 
because of the limitations on the accuracy of computed probabilities. For 
this reason, excessively small values of u and »v (in absolute value) are 
printed, for identification and possible further study of these borderline 
cases. Ifu>0O or v <0, a definite error of at least one unit in the last 
figure of &, is indicated, and the program decreases or increases radii 
by one unit in the last place, according as « > 0,or »y <0, and again 
attempts to bound the correct radius between two numbers of the form 
e~h/2 and¢+h4/2 where c is a rounded number, The checking pro- 
gram provides for the indefinite iteration of this process until the radius 
is bounded in the manner described, but in practice there are only a few 
cases in which the checking program corrects the original &, by one 
unit in the last place, and none in which the correction is more than one 
unit. 


The use of this program requires a detailed knowledge of the 
precision with which ¥ or P is computed, In the high probability 
portions of the tables, since dl’/déK or dP/dk canbe as small as5 x 107 
it was necessary to compute |’ and / to ‘10-% or better, Both tables 
were computed on the NORC, The parameter Jf, which determines 
whether the Taylor or asymptotic series is used, and which is dis- 
cussed in Appendix A, was taken equal to forty. The cpsilon, which 
determines the least upper bound of the smallest term considered, 
whether one series or the other is used, and which is also discussed 
in Appencix A, was set equal to 107'?, 


Program for Checking V(K,c) and P(R,d) Inverse Tables 


ae R, 
rounded value of Ry) 


v@P(a,d)=P 
20 


leeares x 10721 >0, print 2,b,u,Vv; 

If v5 x 1o~tleg, print a,b,u,v; 

If b #0, print a,b,u,v; 

If lul <w, we lul; 

If |v| <w, we fv]; 

Print w at end of entire computation 


h = value of 1 unit in last 

place retained; if 6 dec. digits, 
h = 10°76; R, is rounded by shift- 
ing off extra digits from (Ry+h/2). 
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APPENDIX D 
INVERSE TABLE OF V/K,c) 
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IDENTIFICATIONS FOR THE V(K,c) TABLE 


o,, Sy are standard deviations in 
the x and y directions respectively 
with the mean of the distribution at 
the origin. 


Ko, is the radius of the circle over 
which the integral is evaluated. 
Center of civcle is at origin. 


oe 
¢=— O<se<l (151) 

a. 
4 


This is an inverse table. K is given as a function of V or (P = 
probability) and ec, 


Ranges for the variables. 


Main table: 
c = 0(0.01) 1.00 
(152) 
P = 0,01( 0.01) 0.99 
Supplementary table of high probabilities: 
¢ = 0, ¢ = 0,10(.05) 1.00 
(153) 


P = ,99( ,0005) .9990( .0001) .9999( ,00001) ,99999 ( 000001) .999999 
V(K,e) Example 


Suppose that a group of bombs fall in an uncorrelated bivar- 
jate normal distribution with standard deviations o, = 100 feet, 
o, = 250 feet along the coordinate axes. Ifa point target is at the 
mean point of the distribution, what must the lethal radius of the 
bomb be in order that the probability of destroying the target with a 
single bomb will be P = 0.95? 
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Solution: The roles of x and y must be reversed since it is 
assumed in the theory that Oy § %. Here the parameters are taken as 
o, = 250, Oy = 100, ¢ = oy / Ox = 0.4, P = 0.95, Entering the V(K,e) 
inverse table with these values, it is found that K = R/o, = 2.0051. 
Hence lethal radius R& = 2.0051 o, = 501.3 feet approximately. 
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APPENDIX E 
INVERSE TABLE OF P/R,d) 


NWL REPORT NO. 1768 


IDENTIFICATIONS FOR THE (4, 2) TABLE 


y Distribution is circularnormal 

with common standard devi- 
Ro=R atione, inthe x and x 
directions, and mean at the 
* origin. 

do 
Ro is the radius of the circle 
over which the integral is 
evaluated. Center of circle 
is at (do,o). 


This is aninverse table. R is given as a function of P (probability) 
and d. 


Ranges for the variables 


Main table: 


d = 0(0,1)5( 0,2) 10( 2) 20(5) 120 
(154) 


P = 0,01( 0,01) 0,99 


Supplementary table of high probabilities: 


d = 0, 0.01, 0.10, 0.25, 0.50, 0.75, 1, 1.5, 2, 3, 4, 5, 6, 8, 10, 20, 30, 50, 80, 120 
P = .99(..0005) .9990( .0001) .9999( .00001) .99999( .000001 ) .999999 


(155) 
P(R,d) Example 


Suppose that a group of bomb impact points form a circular 
normal distribution centered at the origin with standard deviation 
0, =0,=0= 300 feet. Ifa point target is at a position given by 
h = 720 feet, & = 960 feet, what must the lethal radius R of the bomb 
be in order that the probability of destroying the target with a single 


bomb be P = 0.60? 


NWI. REPORT NO. 1768 


Solution: Vi? +k? = 1200, d = 1200/0, = 4 


Entering the P(R,d) inverse table with d=40, P = 0.60, itis found 
that R = R/o, - 4.3740. Hence lethal radius R = 4.3740 0, = 1312 feet 
approximately. 
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